Mon. Not. R. Astron. Soc. 000, ??-?? (2012) Printed 20 August 2012 (MN MgX style file v2.2) 



An improved method for the identification of galaxy 
systems: Measuring the gravitational redshift by Dark 
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ABSTRACT 

We introduce a new method for the identification of galaxy systems in redshift sur- 
veys based on the halo model. This method is a modifi ed version of the K-mcans 
identification algorithm developed by lYang et all (|2005h . We have calibrated and 
tested our algorithms usin g mock catalogs generated using the Millennium simula- 
tions (|Springel et al.ll2005f ) and applied them to the NYU-DR7 galaxy catalog (based 
on the SDSS datasets). Using this local sample of groups and clusters of galaxies we 
have measured the effect of gravitational redshift produced by their host dark matter 
haloes. Our results shows radial velocity d ecrements cons i stent with general relativity 
predictions and previous measurements bv lWoitak et al.l (|201l[ ) in cluster of galaxies. 
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1 INTRODUCTION 



In a recent article (jWoitak et all 1201 ll ) (WHH) analyze 
the pattern of spectroscopic redshifts for galaxies in 7800 
clusters from the 7th data release of the Sloan Digital 
Survey (SDSS). Motion of a galaxy in a cluster generates a 
red or blueshift of its spectral lines; the equivalent velocities 
are typically on the order of 600 kms/s. WHH show that 
it is possible to disentangle the superimposed gravitational 
redshift predicted by the General relativity that corre- 
sponds to a radial velocity differences of the order 10 km/s. 
In their analysis the authors used a stacking data from 
several clusters and determined the shift of the redshift 
distribution's centroid with growing radial coordinate. 
The positions and redshifts of cluster galaxies are derived 
from a Ga ussian Mixture B rightest Cluster Galaxy cluster 
catalogue (lHao et al.ll20irjh . one of the largest samples of 
galaxy clusters assembled on the basis of the SDSS-DR7. 
In spite of the fact that massive clusters are expected to 
have the largest gravitational redshift effects, these systems 
are also affected by large g alaxy peculiar velo cities and 
substructure. Kim and Croft (Kim & Croft 2004) suggested 
that it is possible to overcome partially these difficulties by 
averaging over many clusters and groups of relatively low 
mass. 

The existence of volume complete samples of galaxies with 
redshifts measured in the nearby universe as in the SDSS 
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main galaxy sample (|Abazaiian et al] l2009l h lead us to 
explore the gravitational redshift effects using systems of 
galaxies in the nearby universe. Several groups have derived 
catalogs of systems of galaxies applying different techniques 
(Friends of friends FoF, Matched Filters etc.) to large data 
set such as the 2dF and the SDSS. 

In order to make a meaningful comparison between obser- 
vation and theory, group definitions must have a reliable 
three dimensional counterpart. From the point of view 
of current theory of galaxy formation, the most direct 
route for defining galaxy groups is dark matter halos. 
More importan t for this w o rk is t he novel hybrid technique 
introduced by lYang et alj (|2005l ). that starts with cluster 
candidates selected by a classical FoF algorithm and a 
selected sample of bright isolated galaxies. Other satellite 
galaxies are classified according a Maximum Likelihood 
criterion. Such method employs a model of the redshift 
space distribution of the galaxies (as function of system 
mass) as filter. A convolution of the model with the galaxy 
distribution allows to decide the membership of each 
galaxy to the proposed systems (provided by a suitable and 
calibrated FoF method) using an iterative procedure. The 
method relies in an assumed mass to light relation, which is 
a more reliable method for compute the system mass that 
those based on the computation of the velocity dispersion 
and virial r a dii fro m galaxy positions. 

lYang et alj (120051 ) method is a basic version of the well 
known K-means, a clustering algorithm where points are 
assigned to exactly one cluster and all points assigned to a 
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cluster are equals in that system. Yang have carefully tested 
the performance of their group finder. The method was 
applied to the 2dFGRS and compared with those extracted 
from detail e d m ock galaxy redshift surveys. More recently 
lYang et alj i|2007l ) applied this method to the fourth release 
of the SDSS survey in order to study the dependence of 
color, stellar mass, star formation rate and morphology on 
halo mass. 

In order to measure the gravitational redshift caused by 
dark matter haloes in systems of galaxies, we present 
in Section 2 modifications of the algorithm consisting 
in the introduction of a soft degree of assignment of a 
galaxy to each cluster. We have tested and calibrated 
our modified method in Section 3 with an extensive use 
of a galaxy mock catalog building up on the results of 
the Millenium simulations which mimics the SDSS-DR7 
galaxy catalog. In section 4 we use a selected sample of 
clusters and groups of galax ies from the SDSS-NYU-DR7 
datasets l|Blanton et alJ I^OOS*) to compute the gravitational 
redshift produced by its dark matter haloes. This is 
followed by a brief summary. We a dopt for our study a 
flat concordance AC DM cosmology (|Komatsu et al.ll201ll ) 
with Ho = 72km s _1 Mpc,fi M = 0.26, fi A = 0.74, where 
appropiate we define h = if/lOO/km s~ /Mpc. 



2 IDENTIFICATION OF SYSTEMS OF 
GALAXIES: 

The aim of this section is to develop a group finder that 
assigns galaxies in a common halo to a s ingle group. We 
will follow the methodology of lYang et al.1 (|2005T ) which has 
advantages over other group identification algorithms, and 
introduce some improvements related to the K-means im- 
plementation as well as the initial seeds for the model. This 
method relies on the introduction of a background level B, 
that defines a threshold density contrast in redshift space. 
Ideally this value should correspond roughly to the redshift- 
space contrast at the edge of a halo. Taking into account 
theoretical considerations, Yang et al obtain B ~ 10, in ex- 
cellent agreement with the value obtained by maximizing 
the completeness and minimizing the contamination. 
Nevertheless, we first notice that it is not necessary to intro- 
duce such a high threshold galaxy density contrast. Rather 
we allow a lower background level value so that the likeli- 
hood of a galaxy in a given halo is greater than that of the 
main or central galaxy of another nearby haloes. Such value 
was chosen using the simulated catalogs described below. We 
can therefore use the likelihood value as main/satellite clas- 
sification criteria by classifying a galaxy as satellite if the 
likelihood of pertenence to another halo is over the back- 
ground level. We argue that this methodology could be used 
as a classification process. Rather than starting from pre- 
fixed seeds halos (such as in the work of Yang et al.), each 
galaxy in the catalog is allowed to start as central galaxy in 
the halo where it resides. As the iterative method proceed, 
many galaxies subsequently will be classified as satellites. 
No further criteria for the seeds selection are needed with 
this improvement. This avoid the problems introduced by 
th e pre-selection of see ds using a FoF algorithm ( discussed 
bv lBerlind etafl (|2006l )). 



2.1 From hard to soft assignments: 

Yang's cluster assignment method can be considered as an 
example of a competitive learning algorithm: the well known 
K-means. The K-means is an algorithm for assigning N data 
points (i.e. galaxies) into K clusters (i.e. groups of galaxies). 
Each k cluster is parameterized by a vector Q k (which is 
called the cluster mean) with elements position of the cen- 
ter and size of the group of galaxies. The algorithm works by 
allowing the K-clusters to compete each other for the right 
to own the data points. This procedure owes its advantage 
in its simplicity, and in its local-minimum convergence prop- 
erties. However it as a main shortcoming, is a "hard" rather 
than "soft" algorithm: points assigned to a cluster are equals 
in that cluster. Points located near the border between two 
or more clusters should also play a role in determining the 
locations of all the clusters that they could plausibly be as- 
signed to. Nevertheless in the K-means algorithm, each bor- 
derline point is dumped in one cluster, with the same vote 
than all the others points in that cluster, and no vote in 
any other cluster. The previous criticisms of K-means mo- 
tivate th e "Sof t K-means algorithm" used in this work (see 
iMacKavl (|2003ft ). 

A fundamental hypothesis for this clustering problem is the 
underlying idea that the maximum likelihood method, which 
identifies the setting of the parameter vector H that max- 
imizes the likelihood P(Data|fi, Model) provides a good fit 
to the data. Briefly, the algorithm is set in three steps: Ini- 
tialization, Assignment and Update. 

Step 1: Initialization In order to start the K-means, the 
mean values Q.' k:i must be initialized in some way, we start 
using all the galaxies in the catalog. A representation of this 
assignment of galaxies to groups is given by the so called 
"responsibilities", which are indicator variables that a 
given galaxy n belongs partially to the k group of galaxies. 
We firstly take all galaxies as 'central', and so are considered 
as the center of a potential dark matter halo. Using the total 
group luminosity L groU p and a model for the group mass- 
to-light ratio we can estimate the halo mass associated i.e. 
the halo radius r2oo, the virial radius r v i r , the virial velocity 
Vvii = (GM/rvh ) 1 / 2 and the line-of-sight velocity dispersion 
of the galaxies within the dark matter halo a = VW/v2. 
The total luminosity of a selected potential group is esti- 
mated using (at this stage is simply the luminosity of the 
galaxy). L gIoup — y^. r^Lj, where L; is the luminosity of 
each galaxy in the group, r£*^ is the k group responsibility 
for each (i) galaxy (=1 initially). 

Step 2: Assignment step Once we have a group cen- 
tre, and a tentative estimate of the group size, mass, and 
velocity dispersion, we can assign satellite galaxies to this 
group according to the properties of the associated halos. 
If we assume that the phase-space distribution of galaxies 
follows that of the dark matter particles, the number den- 
sity contrast of galaxies in redshift space P(Data|fi fc , Model) 
around the k group centre at redshift 2 gr0 up is computed in 
a similar way than in Yang's technique. 

P(Data|ff fe , Model) = P M (R,Az) = ^H^lp(Az) , (1) 

c p 

Here Az = z — % r ou P and £(-R) is the projected sur- 
face density of a (spherical) NFW halo. Hereafter R is the 
projected distance at the redshift of the group. The func- 
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tion p(Az)dAz describes the redshift distribution of galax- 
ies within the halo for which standard Gaussian shape is 
adopted, with a as the rest-frame galaxy velocity disper- 
sion. Thus defined, Pm (R, Az) is the three-dimensional den- 
sity contrast in redshift space. In order to decide whether a 
given galaxy should be assigned as a primary or as a satel- 
lite into a particular group, for each galaxy we loop over 
all groups, and compute the corresponding distance (R, Az) 
in phase space between galaxy and group centre. Firstly we 
classify the galaxies as primaries or satellites according to 
the value of P(DatajH, Model) starting from the most mas- 
sive systems towards the smaller ones, as was described in 
subsection 2.1. Instead of assigning each galaxy (n) to the 
group for which Pm{R, Az) has the highest value, we pro- 
ceed to evaluate the responsibilities for each group (k) as: 



(n) 7r (fc) P(Data|f2 fc , Model) 
J2k 7r (£) P(Data|n*, Model) 



(2) 



Step 3: Update step In the update step, the following 
model parameters (" ): center angular position and group 
redshift, are adjusted to match the sample mean of the data 

points x (n> that they are responsible for = — 



where = rj: is the total responsibility of mean 

k. Now the algorithm also includes groups weighting pa- 

r(fc) 



rameters 7Trfci = 



r which also update themselves, al- 

lowing accurate modeling of data from groups of unequal 
weights. The total responsibility of a cluster becomes a bet- 
ter measure of the real occupation of the halo, and being a 
continuous parameter it is more adequate for modeling and 
interpretation. Since some of the groups parameters (mass, 
velocity dispersion, virial radii) are derived from the total 
group luminosity, we could made use of the responsibilities 
in order to obtain a "soft" weighted luminosity, which should 
be more representative of the real value. This procedure has 
the advantage that the inferred luminosities (i.e. the masses) 
of the groups are almost independent of the possible con- 
tamination by interlopers. It should be recalled that the 
memberships of the selected groups are remarkably insen- 
siti ve to the adopted mass-luminosity model, as was stated 
by (|Yang et alj|2005l ). Next, we iterate between Step 2 and 
3 until there are no further changes in the memberships. 
Finally we assign membership for each galaxy to the group 
which have the largest responsability (in case of a tie, we 
assign to the most massive group). 

Note that, unlike with the usual FOF method, this group 
finder also identifies groups with only one member. Since we 
start using all the galaxies in the catalog, we could poten- 
tially assign each galaxy to their corresponding dark matter 
halo, since we include all possible groups seeds. 



out so far. In what follows, we will use the semi- 
analytic galaxy catalogues at redsh ift zero cons t ructed 
from the Millennium simulation by ICroton et alj (|2006T ) 
(http:/ /www. mpa-garching.mpg.de/galform/agnpaper/ ), 
who considered a detailed model for cooling, star formation, 
supernova feedback, galaxy mergers and metal enrichment 
as well as a simple treatment of heating by a central AGN. 
These authors find that several observed properties of 
the galaxy population at z — are suitably reproduced 
indicating that galaxy characteristics and stellar ages are 
in much better agreement with observation than for models 
without AGN feedback. We build up mock catalogs of the 
SDSS-DR7 (spectroscopic) survey in order to calibrate 
our system identification method as well as comparing 
observed galaxy properties of models and observations. 
Since the median redshift of galaxies in SDSS is ~ 0.1, the 
r-band absolute magnitudes M r of each model galaxy are 
K-corrected to z = 0.1 (M®' 1 ) using the K— corre ction code 
(kcorrect v3_lb) of iBlanton fc Roweis! {2007) and the 
luminosity evolution model of Blanton et al. (2003). Galaxy 
redshifts are assigned by placing the observer at the corner 
of the simulation box and are determined by the comoving 
distance to the observer plus the galaxy peculiar velocity 
contribution. The corrected r- band magnitudes are given 
by: Mr' 1 = — 2.5xlogL + K corre ction + E correc u n — 51og/i. 
The survey geometry and the radial selection function 
were properly taken into account in the construction of the 
mock catalogue where we record galaxy redshifts, angular 
coordinates, 3D positions, apparent and absolute u,g,r,i and 
z band SDSS magnitudes, as well as an index to identify the 
galaxy assignment to a dark matter halo. In order to model 
the mass-weighted luminosity relation we select a common 
luminosity scale L19, defined as the luminosity of all group 
members brighter than M®' 1 = — l9 + 5*log(h). To calibrate 
this relation for the weighted luminosities we select from the 
mock SDSS catalog all groups with z < 0.068, correspond- 
ing to the volume limited magnitude given for a catalog 
with apparent magnitude limit of 17.5. Using the actual 
group centers and masses we compute the responsibilities 
of each group over all the galaxies in the catalogue. The use 
of responsibility weighted luminosities provides a slightly 
improved (with lower disp e rsion) mass-luminosity relation, 
compared to lYang et al.l l|2005T ). Since the total group 
luminosity is dominated by the few brightest galaxies, the 
estimated mass is much less sensitive to the absence of faint 
group members (typical in flux-limited catalogs) and the 
effects of interlopers. This methodology allows to avoid one 
of the problem of the maximum likelihood, the well known 
phenomenon of overfitting. For all groups selected below 
the redshift limit 2 < 0.068 we compute L19 directly from 
the selected members with M r )A — 5 * log(h) —19.0. For 
groups at higher redhifts, we compute L SIOup and use the 
average relation between Lig and L SIOup to estimate the 
former. 



3 APPLICATION TO SIMULATED AND SDSS 
CATALOGS 

Numerous mock catalogs have been produced from 
full-bl own semi-analytic model of galaxy forma- 
tion (Wang ct al. 2006). The Millennium Simulation 
jSpringel et all l2005h used in this work, is one of the 
largest simulation of cosmic structure growth carried 



3.1 Performance of the method 

A real group is defined as the set of galaxies that reside in 
the same halo of the mock catalog. In order to quantify the 
group finder performance we introduce the completeness, 
defined as the ratio between the number of true members 
selected by the group finder and the total number of true 



4 Mariano Javier de Leon Dominguez Romero 



i o ° i o 1 10 J io 3 io 10 ' 

Viria Moss [ h _1 72 M ] 



10 I0 10"10 I2 10 I5 10 I1 10 15 
Virial Mass [ h~' « M.] 



Figure 1. Left panel: completeness (open triangles) and contami- 
nation (open squares), Right panel: fragmentation (open squares) 
and purity (open triangles) of groups (averages) as a function of 
halo mass for systems identified in the mock catalogue based on 
the Millennium simulations. Errorbars correspond to 1-cr vari- 
ance. 



group members, and the contamination, defined as the ratio 
between the number of false members (interlopers) selected 
by the group finder and the total number of real members. 
We also measure the fragmentation, defined as the mean 
number of extra galaxy groups per dark matter halo having 
mass at least 0.1 times that of their true associated galaxy 
group and compute the purity of the systems defined as the 
ratio between the true members and the group members. 
Figures Q] gives an overview of the completeness, contami- 
nation, fragmentation and purity for our group finder. Dots 
correspond to individual groups (halos), while large triangles 
indicate the average for all the systems included inside the 
window function without restrictions in the absolute magni- 
tude limits. As can be observed, for high mass/population 
systems it is difficult to recover all the members in the 
outer regions. Nevertheless, these objects are preferentially 
infalling galaxies that are difficult to recover by any method. 
Fragmentation of groups should be considered if these sub- 
groups comprise a considerable mass fraction. If this frac- 
tion is ~ 10% the effect is negligible for high mass systems 
but affects intermediate mass/occupation systems. With the 
modifications introduced, the improved algorithm is signif- 
icantly more efficient than the popular FoF identification 
algorithms or its variants in terms of high completeness and 
purity, low contamination and fragmentation as is shown in 
Figure [T] Notice that our test are based on a mock sam- 
ple from a full N-body SAM based model, which presum- 
ably results in a more realistic scenario than a CLF base d 
one as is the 2dF mock catalog used by I Yang et all |2007l ). 
Three important output parameters of the method should be 
stressed: each group responsibility for a galaxy provides an 
effective medium for interlopers classification (as an exam- 
ple a galaxy with similar responsibilities by two different and 
nearby haloes); the total responsibility of each halo, which 
could serve as a new halo occupation measure (introducing 
some limit on magnitude); and the total weighted luminos- 
ity as an important intermediate step to the computation of 
the total dark matter halo mass. 



3.2 System Identification in the SDSS-NYU-DR7 

Using a dedicated 2.5 meter telescope, SDSS provide the 
largest and most complete photometric and spectroscopic 
galaxy survey available. We have used the large-scale struc- 
ture sample samplel4 fr om the NYU Value Added Galaxy 
Catalog (NYU-VAGC; iBlanton et all (|2005l )) DR7 version 
as our primary galaxy sample, a complete sample of galax- 
ies with reddening corrected r magnitudes brighter than 
17.5. We take into account fiber collisions by giving each 
collided galaxy their photometric redshift estimate. In order 
to limit the effects of incompleteness on our group identifica- 
tion procedure, we restrict our sample to regions of the sky 
where the completeness (ratio of obtained redshifts to spec- 
troscopic targets) is greater than 90%. Galaxy magnitudes 
are corrected for Gal actic extinction using the dust maps of 
( Schlcgcl et al . 1998), absolute magnitudes are k-corrected 
|Blanton fe Roweisll2007r> and corrected for passive evolution 
jf31anton et alJ l2003h " to rest-frame magnitudes at redshift 
z — 0.1. We use the fit Q(z) = 1.6(2 — 0.1), which takes into 
account the effects of the luminosity evolution. We select 
all galaxies with extinction corrected apparent magnitude 
brighter than r = 17.5, redshifts in the range 0.015 ^ z ^ 
0.15, and with redshift measurement completeness ^ 0.9. 
We use a lower redshift limit z > 0.015 to alleviate some of 
the problems asociated with the accuracy of the photometry 
of some nearby extended galaxies. For a uniform magnitude 
limit, our selection function n(r), can be separated into the 
product of an angular and a radial part: n(r) = n(r)n(r), 
where R = rf and f is a unit vector. The angular part may 
take any value between and 1, and provides the complete- 
ness as a function of position, i.e., the fraction of all sur- 
vey selected galaxies for which reliable redshifts are given. 
In order to compute the selection function, we use Mangle, 
http: //casa.colorado.edu / ajsh/mangle. Such codes 
ijHamilton fc Tegmarkl |2004| ). are designed to deal accu- 
rately and efficiently with co mplex angular mask. We ap- 
ply the procedure suggested bv lBerlind et al.l (|2006l )) to the 
groups identified by our algorithm on the NYU-SDSS-DR7 
galaxy sample using the mask in order to avoid completeness 
issues due to the angular footprint. 



4 A MEASURE OF GRAVITATIONAL 
REDSHIFTS 

Using the systems of galaxies identified in the SDSS-NYU- 
DR7 galaxy redshift survey we measured the gravitational 
redshift due to the dark matter gravitational profile. A stan- 
dard stacking technique is used in order to disentangle the 
kinematic Doppler effect from the gravitational redshift, 
given that the latter shifts the centroid of the observed ve- 
locity distribution. Coordinates and redshifts of cluster cen- 
tres need no corrections since our identification methodology 
provides reliable values. We search for all galaxies within a 
window in phase space of 3 virial radii and 3cr in velocity 
around the cluster centres. As the final step we combine 
redshift data of all clusters into one. Since we are interested 
in the detection of the GR effect in intermediate mass sys- 
tems, we have considered two samples of systems of galaxies: 
low mass (1O 13 M to 1O 14 M , groups of galaxies) and high 
mass ( 1O 14 M0, comparable to the cluster sample analyzed 
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by WHH). Our final sample comprises systems of galaxies 
in the range 0.015 ^ z ^ 0.15 with a mean redshift of 0.1 
hosting on average 28 galaxies with spectroscopic redshift 
measurements for clusters and 14 for groups. Those groups 
with less than 4 redshift determinations were not included 
in the sample. Since the presence of interlopers in our sam- 
ple is negligibly small, it is possible to measure directly the 
mean value of the velocity distribution of cluster galaxies, 
A. Assuming spherical symmetry and no strong inhomo- 
geneities, the gravitational redshift profile of a system can 
be calculated using the following formula i|Cappilll995h 



A S (R) 



c£(R) 



Mr) ~ *(0)] 



p(r)rdr 



R 2 



(3) 



where R is the projected cluster-centric distance, <E>(r) is the 
gravitational potential, p(r) and are the 3D and sur- 

face (2D) density profiles of galaxies. In order to estimate 
this effect for the data combined from a cluster sample, it is 
needed to convolve this expression with the distribution of 
cluster masses in the sample. In order to test the possibility 
of determining this effect in our sample we introduce the cor- 
responding blueshift in the mock catalog of galaxies in their 
parent haloes. Given that the dark matte r haloes follow an 
NFW density profile l|Navarro et al.|[l996l ). the gravitational 
potential results: 



*(r) = -(GAfv/rv) 1 '^) 1 ' 2 ^ + r/rv) 

r/r v 



(4) 



and g(c v ) = l/[ln(l + c v ) — c v / (1 + c v )] with c v the concen- 
tration depending on the halo mass. Using formulas 3 and 4 
we introduce the corresponding gravitational redshift effect 
by adding A to the line of sigth velocities of the member 
galaxies (Hubble flow plus peculiar velocities) that popu- 
lated each dark matter halo in the Millennium simulation. 
In the left panel of figure [5] we show A as a function of the 
projected radial distance and its comparison to the imposed 
effect in the systems identified in our mock catalogs. This 
should be compared with the measurement showed in the 
right panel for the SDSS-NYU-DR7 systems. Upper/Lower 
panels correspond to low/high mass systems. Although the 
predicted signal is significantly smaller for groups than for 
clusters of galaxies (showed in dashed lines), the quality and 
size of the sample of low mass systems made possible to 
clearly detect the gravitational redshift effect predicted by 
the General Relativity. We notice that the high mass results 
are in agreement with jWoitak et al.ll201ll ). 



5 CONCLUSIONS 

• An efficient halo finding algorithm has been developed 
and tested using mock catalogs derived from numerical sim- 
ulations with semianalytic galaxy formation methods. The 
performance of the method is solid and efficiently provides 
complete, almost bias free galaxy system catalogs based on 
galaxy redshift surveys. 

• The new method was used to compile samples of groups 
and clusters of galaxies in the SDSS-NYU-DR7. The clusters 
of galaxies sample was used to reproduce the results found 
by WHH in clusters of galaxies. Following the suggestion by 
jKim fc Croftll2004h . we are able to detect the gravitational 
redshift in the sample of groups of galaxies. 



[i 
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Figure 2. Observed gravitational redshift (points with 1-cr error- 
bars)in terms of the mean velocity A, as a function of the pro- 
jected clustercentric distance for Millennium simulated systems 
in the mock catalog (left panels) and SDSS-NYYU-DR7 systems 
of galaxies (right panel). Upper and lower panels correspond to 
groups and clusters of galaxies respectively. Dashed lines show 
the General Relativity theoretical predictions. 



• Our measurements of the gravitational redshift effect 
by dark matter haloes show a good agreement with the pre- 
dictions of General Relativity. 
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